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ABSTRACT 

A five-level baroclinic model for numerical prediction of height 
fields at 1000, 850, 700, 500, and 300 mb is developed and programmed. 
Several forms of the model are tested to determine the contribution 
of the various terms of the vorticity equation. Forecasts from three 
different methods are compared. The model shows utility in prognosis. 
Although there are individual differences among the methods, no one of 
the three is significantly better over the entire hemisphere on the 
data tested here. 

The writers wish to express their appreciation to Professor George 
Haltiner of the U. S. Naval Postgraduate School for his assistance and 
guidance in this investigation. Appreciation is also expressed to the 
personnel of the U. S. Navy Fleet Numerical Weather Facility for their 
cooperation during the preparation of this paper. We are especially 
indebted to CDR L. C. Clarke and LT G. Lawniczak for use of their 
solution of the omega equation, and for their assistance to us, and 
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1. Introduction. 

Until very recently, the barotropic models have been the principal 
source of numerical weather prognoses. Their shortcomings in areas of 
baroclinic development, plus the advent of sophisticated digital computer 
systems, have spurred interest in baroclinic models. [It is the purpose 
of this research to develop and program a five-level baroclinic model and 
to investigate the contribution of the various terms of the vorticity 


equation to the prognosis of the pressure field. 





Die Background. 
In the model developed here, prognosis is accomplished by means of 
the vorticity equation, which, with pressure as the vertical coordinate, 


can be written 
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Basic to the derivation of this equation is the Helmholtz theorem which 
expresses the wind as the sum of a divergent and a non-divergent compon- 


ent: 
Velk x VW + \/X, (2) 


The potential function X for the divergent part of the wind is obtained 


through solution of the continuity equation 


fe ° 
ae oa ae” (3) 


The stream function W for the non-divergent portion is approximated by 


solution of the balance equation 
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as described by G. ‘Arnason eal where Od is the geopotential. Use of 
the stream function in the vorticity equation enables one to avoid the 
geostrophic approximation and the concomitant spurious anticyclogenesis. 


The vertical component of relative vorticity, ¢ , is the Laplacian 


of the stream function 
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and absolute vorticity, // , 1s the sum of this and the vertical component 


erthe earth si vorticity.,.f. 


Yr (6) 


Vertical velocity,CO, is computed by the diagnostic equation 
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developed by Clarke and Lawniczak [3]. The forcing function is calculated 
with the geostrophic approximation. This results in some inconsistencies 
in the model but to avoid reprogramming the entire equation, we have gone 
ahead with comparison-type experiments. Vertical velocity is assumed zero 
at 100 mb, but other boundary conditions are among the options available 
to the user of this computer program as follows: 

a) the lower boundary can be level or include terrain; 

b) vertical velocity can be set equal to zero at the lower boundary, 
or computed with or without terrain effects and with or without the 
effect of friction; 

c) the stability parameter ¢ can be a point variable, a function of 
pressure alone, or a constant; and 

d) the maximum acceptable residual in the 3-dimensional relaxation 
can be selected by the user. 

In the vorticity equation, the contributions from the twisting term, 


O. o> p° 
(22 =r — 22 as and the vertical advection of vorticity term, 
J 
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Pp? are usually considered small, and some investigations have 
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shown that they tend to compensate each other GE Also, integral theo- 


rems show [4] that to avoid the mean generation of vorticity in a closed 





volume, the two terms should be included or neglected together. Both 
have been omitted from the initial investigation reported here. 
A further simplification involves omission of the advection of abso- 


lute vorticity by the divergent portion of the wind, 2% alt Sis ee ay) ; 





For consistency, it is necessary then also to substitute the mean value 
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of the Coriolis parameter, { , for the absolute vorticity in the diver- 
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gence term, 


In the equations that follow, the stream function has been redefined 


so as to be in units of height, 


Y= | (8) 


The resulting forms of the vorticity equation, expressed in finite- 


difference notation, are 
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and the complete equation 
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where Ci is the grid distance and /W is the ratio of this distance to the 


corresponding geographical distance. 





Horizontal spatial operators are approximated by the usual centered 
differences. To approximate the vertical derivatives in the vorticity 
equation a Lagrangian parabola was fitted through three points. For a 


centered difference the equation is 


DA. Cet (p-PP) y, CPP, 
: x oa ) A. panels! ae (12) 
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All computations were done on the Control Data Corporation model 
1604 computer, using the U. S. Navy Fleet Numerical Weather Facility 
octagonal grid. This grid, centered at the pole and inscribed in the 10N 
latitude circle, consists of 1977 grid points with a grid size of 381 km 
at 60N. 

The week of 22 November 1962 was a period of considerable cyclogene- 
sis in dense-data areas, and therefore excellent for testing a baroclinic 
model. Initial data used were the U. S. Navy Fleet Numerical Weather 
Facility height analyses at 1000, 850, 700, 500, and 300 mb for 12002 on 
22 November 1962. The height analyses for 0000Z and 1200Z 23 November 


were used in the evaluation of the predictions. 


Dis Procedures. 

In the diagnostic G@-equation for the initial tests discussed here, 
the boundary conditions used are: 

a) vertical velocity equal to zero at 100 mb; 

b) vertical velocity equal to zero at the lower boundary; 

c) no terrain effects at the lower boundary; and 

d) the stability parameter, , , is a function of pressure. The pro- 
gram is used in "forecast mode''; that is, temperatures are computed hydro- 
statically, rather than taken from objective analysis. The maximum 
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residual in the 3-dimensional relaxation is 1 x 10 “ mb per second. 

The sequence of calculations programmed can best be explained with 
reference to figure 1. Input data in every case are the U. S. Navy Fleet 
Numerical Weather Facility height analyses at the five standard levels, 
modified to ensure vertical consistency. 

In the first method with equation (9), stream functions are used 
throughout. The inversion to height fields [1], effectively the reverse 
of the balance equation, is thus bypassed and stream functions are used 
for input data for both diagnostic and prognostic equation. This is 
hereafter referred to as method I. The resulting 24-hour prognoses are 
displayed on charts 6 through 10 and discussed in section 4 below. 

Since the geostrophic approximation is used in the forcing function 
of the W-equation, and temperatures are computed hydrostatically, there 
might be some reason to suspect the vertical velocities computed from 
height fields could be more accurate than those from stream functions. A 
completely geostrophic model was programmed, using pressure-height data 
throughout. This method showed the expected spurious anticyclogenesis 
even at 12 hours and was carried no further. 
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The complete process diagrammed in figure 1 has been used in the re- 
mainder of this investigation. In the first time step the initial data 
(five vertically-consistent height fields) are used as input to the bal- 
ance equation and to the &)-equation. From the former are obtained stream 
functions for the five levels; the latter yields vertical velocities at 
the same levels. With these stream functions and vertical velocities the 
vorticity equation (equation (9) for method II, equation (10) for method 


LE 
i 
III) gives re 





, the local change in stream function, at the five levels. 
When these are added to the initial values, the first forecast time step 
is complete. As shown in figure 1, these new E values are now input 
fields for both the vorticity equation and the balance equation; the 
latter being used for inversion to pressure height as required by the 
() -equation. 

Subsequent time steps follow the same routine, except in the applica- 


tion of SE. which is done according to the centered-difference equation 


Gapopebe eat Bia (13) 
In order to minimize the errors inherent in the initial forward 

time step, the first time increment is only 15 minutes. Centered time 

steps are then used to complete successive 15-minute prognoses to one 

hour, and finally l-hour prognoses to 12 hours. The fields are then 

smoothed lightly to remove any small perturbations resulting from trunca- 


tion or round-off errors, and the process is repeated to 24 hours. 





4. Results. 

In review: method [ uses stream functions throughout, and with 
equation (9) as the prognostic equation excludes the advection of abso- 
lute vorticity by the divergent portion of the wind and substitutes ic 
for » in the divergence term; method II inverts the stream-function 
fields to pressure heights at every time step, but retains the other 
simplifications of method I; method III also inverts to pressure height 
at each time step, but uses the more complete equation (10) as the prog- 
nostic equation. 

Charts 1-5 are the pressure height analyses at the beginning of the 
period. The next 15 charts are the results of the three 24-hour fore- 
casts: method I in charts 6-10; method II, charts 11-15; and method III, 
charts 16-20. Charts 21-25 are the observed height analyses at the end 
of the forecast period. Chart 26 is the U. S. Navy Fleet Numerical 
Weather Facility barotropic forecast at 500 mb for the same period. 
Table 1 shows the contour origin and interval for all the charts. 

In general, corresponding charts from the three methods are quite 
similar. They all retain a typical meteorological appearance and the 
various systems show a reasonable vertical structure. Further indica- 
tion of their similarity is contained in table 2. The "pillow" referred 


to is 


\A77 


_2, (A- Ee 


Pillow =) \-W—s (16) 
1477 


and the RMSE is the root-mean-square error calculated in the usual 


fashion: 





; (17) 





Nevertheless, a level-by-level comparison of the three forecast 
methods shows some differences among them. At 1000 mb, for instance, the 
two low centers over North America merged and deepened, while the trough 
moved nearly 20° east. All methods forecast the merging and consider- 
able trough movement, with method [I giving the best placement of the 
trough. During the forecast period a low center developed in the area 
of weak gradient southwest of the Kamchatka peninsula. This development 
was forecast by all methods, although with too much cyclogenesis. 

Method I was again the best. In no case is there any evidence of spuri- 
ous anticyclogenesis, even at low latitudes. 

The forecasts at 850 and 700 mb show movement consistent with that 
forecast at 1000 mb. By 700 mb the tilt of the baroclinic systems has 
become apparent. A good example is the stationary low center over the 
Northwest Territories, which deepened by 26 meters. All forecasts held 
it stationary, but methods I and III forecast two and seven meters deep- 
ening, respectively, while method II forecast a 2-meter filling. The 
weakening of the low center east of Greenland was forecast by all methods, 
but best by method II. 

At 500 mb, since the chart of the U. S. Navy Fleet Numerical Weather 
Facility barotropic forecast has been included, there are four charts to 
compare. There is a striking similarity between the barotropic progno- 
sis and that from method I. The method II and III prognoses, themselves 
much alike, are quite different from the other two. Method [I and the 
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USNFNWF model both forecast the movement and reorientation of the trough 
over Eastern United States, and the weakening and displacement of the 
low center near Hudson Bay. Methods II and III show the trough lagging 
far behind its true position, and retain the low as a stationary system. 
Method III even deepened it slightly. In contrast, methods II and III 
produced much better forecasts than the others in the western Europe- 
Mediterranean area. 

At 300 mb all three methods agree in general on the movement and 
development of the major systems. This development is in some cases out 
of phase with the actual development, but no one method can be selected 
as better than the rest. The initial and verifying charts at 300 mb 
show more effect from the process of making the data vertically consist- 
ent than was evident at lower levels. 

Referring again to table 2, there seems to be no best method among 
the three tested. Certainly the addition of the advection of absolute 
vorticity by the divergent wind and the use of absolute vorticity in the 
divergence term produced little change and no over-all improvement in 
this particular forecast. The effect of the inclusion of the vertical 
advection of absolute vorticity and the twisting term remains to be 


tested. 
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Table 1 


Level Contour origin Contour interval 
(mb ) (meters) (meters) 
1000 120 60 

850 1440 60 

700 3060 60 

500 5580 60 

300 9120 120 


12 





Table 2 








Method Level Pillow RMSE 
(mb ) (meters) (meters) 

I 1000 ~ey.. 9 43.6 

850 - 1.8 37.6 

700 - 1.8 41.8 

500 - 1.2 DES 3 

300 13.4 80.3 

II 1000 Lee 43.3 

850 SZ at 40.8 

700 0 42.7 

500 ao 53.4 

300 ee, 1d 

III 1000 Seo 49.1 

850 ) - 3.4 43.0 

700 eet 43.6 

500 sc) 54.4 

300 13.4 80.0 

USNFNWF 500 - 4,9 47.0 

Persistence 500 =353 57.4 
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APPENDIX I 


The following charts were computed during the course of this study. 


They are included here for easy reference. 


Charts 1-5 Analyses, 12Z 22 Nov 1962 
Charts 6-10 24-hour prognoses, method [I 
Charts 11-15 24-hour prognoses, method II 
Charts 16-20 24-hour prognoses, method III 
Charts 21-25 Analyses, 12Z 23 Nov 1962 
Chart 26 24-hour prognosis, USNFNWF 
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